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In a superconducting ferromagnet, the superconducting state appears in the ferromagnetic phase 
where usually a domain structure has already developed. We study the influence of the super- 
conducting screening currents on a stripe structure with out-of-plane magnetization, in a film of 
arbitrary thickness. We find that superconductivity always induces a shrinkage of the domains, 
and there is a critical value of penetration depth below which a mono-domain structure is more 
stable than the periodic one. Furthermore we investigate the possible different effects of singlet 
and triplet superconductivity on the domain width, as well as the conditions for the existence of 
vortices in the domains. The obtained results are then discussed in light of the experimental data 
of superconducting ferromagnets URhGe, UGe2, and UCoGe. 
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I. INTRODUCTION AND MODEL 

In the discovered Uranium-based superconducting fer- 
romagnets (SFM) UGeA URhGe 2 . and UCoGe^ the 
Curie temperature 9 is much higher than the supercon- 
ducting critical temperature T c which means that super- 
conductivity appears in the ferromagnetic state where it 
is usually a domain structure (DS) which develops. In 
previous works^r— the influence of the superconducting 
screening currents on the DS has been studied in the 
case of thick or bulk systems, when the thickness 2L Z 
along the easy magnetization direction is much larger 
than both the transverse domain width I (see Fig.[TJ and 
the London penetration depth A. The domain period 21 
at equilibrium results from the balance between a pos- 
itive contribution to the energy density due to domain 
walls and a negative contribution from the magnetic in- 
duction, as the magnitude of both increases when I is 
reduced. An exact energy minimization 7 shows that for 
the condition A > w/(8ir) (where w is an effective domain 
wall width that parametrizes the wall energy) supercon- 
ductivity decreases the domain size due to partial pene- 
tration of the magnetic field near the domain wall. This 
energy decrease is proportional to A and the formation of 
the domain wall is favorable when this contribution coun- 
terbalances the energy of the domain wall itself, which is 
proportional to w. For A < u?/(87r) the system is in a 
mono-domain state without any domain wall^— . In the 
present article we extend Faure and Buzdin's work 7 - by 
deriving the expression of the energy valid for any thick- 
ness 2L Z and we discuss the DS in all limits, in particular 
when A, I 3> L z . Furthermore we investigate the possible 
effects of singlet and triplet superconductivity on the DS, 
as well as the conditions for the existence of vortices in 
the domains. 

Model. — We consider a ferromagnetic film of thickness 
2L Z that can also become superconducting (see Fig. [TJ. 
The z-axis is chosen perpendicular to the film with the 
surface edges at z = ±L Z . Domain walls parallel to the 
j/z-plane separate the periodic ferromagnetic structure 
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FIG. 1. (color online). Geometry of the considered stripe 
domain structure: the domain width is I and the film thickness 
is 2L Z . The alternating magnetizations M of the domains are 
perpendicular to the film. 



into domains of equal width I and with magnetization 
M = M(x)e z alternating along the x-axis, i.e. M(x) — 
±M Q = (4M /0 Er=o sm{qx)/q with 
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This means that we consider the domain wall thickness 
very small compared to I and A. The energy per surface 
unit is J"(B, 1) = J" M (B, I) + .Fsc(B, I) + Jbw(0 where 
the magnetic field energy is given by the relation 
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and the superconducting current energy is expressed in 
the London limit 9 (i.e. the superconducting coherence 
length £ < A, I) 
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Here, d x d y is the total area the film surface. The explicit 
dependences with respect to the domain width result 



from the Ansatz we made for the magnetization M(Z). 
The contribution from the domain walls may be written 
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where w is a domain wall characteristical length scale. 
Usually w provides an upper limit of the real domain wall 
width 7 . Hereafter, the temperature T will be considered 
phenomenologically within an explicit dependence of the 
London penetration depth X(T). The latter is finite in 
the superconducting low— T phase, and diverges at the 
critical temperature T c . 

Outline of the article. — In the second section we ana- 
lytically calculate the field distribution and the energy of 
a periodic stripe DS in a SFM of arbitrary thickness. An 
exact expression of the energy is obtained as an infinite 
sum. In section III the equilibrium size of the domains 
is determined in the superconducting state as a function 
of the penetration depth. Analytical expressions of the 
domain width are obtained in limit cases. The stabil- 
ity of the domain structure is also discussed. We then 
investigate the difference between singlet and triplet su- 
perconducting states. Specifically, the triplet pairing is 
described by a domain wall energy which is temperature 
dependent, in contrast with the singlet pairing. In sec- 
tion IV we draw the condition for vortices appearance at 
the center of a domain. Finally in section V we apply 
our results to the DS in Uranium-based superconducting 
ferromagnets before concluding. 



II. METHOD OF SOLUTION 

First, we minimize the energy J-"(B,Z) with respect to 
the magnetic induction B. This yields the London equa- 
tion A(B - 4ttM) = A~ 2 B in the film and AB = 
outside. Once B(Z) is found by solving the London equa- 
tion, the resulting total energy, J"(Z) = J r (B(Z),Z), will 
be minimized with respect to I to determine the domain 
width at equilibrium. 



A. Magnetic field B(7) 

We find the magnetic induction B by Fourier expan- 
sion. Using Maxwell-Thomson equation V • B = and 
symmetry relations B z (—x) = —B z {x) and B x (—x) — 
B x (x), the field can be written as 
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where we remind that q = (2k + 1)tt/1. When solving 
the London equation in the film and Maxwell's equation 
outside, we use the symmetry condition B z (—z) = B z (z) 
and the continuity condition at the surfaces. This yields, 




FIG. 2. (color online). Magnetic field distribution of a wide 
domain (I 2> L z ) (a) in the normal state and (b) in the super- 
conducting state (A = 3L Z ). Flux lines (gray-and- white solid 
lines) and the contour plot of the field magnitude |B|/47rMo 
(shaded background) are shown between x = and x — 1/2. 
The horizontal dotted line is the upper film surface. 



for \z\ < L z , 

A ( ) = JL (l qcosh{q z z) 
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and for \z\ > L z 
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Figures [2] and |3] show the distribution of the magnetic 
field in the normal state and the superconducting state. 
It is plotted between the domain wall and the domain 
center. In the normal state, the magnetic field distribu- 
tion in a wide domain (I 3> L z ) is concentrated around 
the domain wall (see Fig.[2][a)) while in a narrow domain 
(/ <C L z ) it is nearly uniform and equal to ±4irMoe z (see 
Fig-GUa)). The supercurrent screens the field around the 
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FIG. 3. (color online). Magnetic field distribution of a nar- 
row domain (I <C L z ) (a) in the normal state and (b) in the 
superconducting state (A = 0.02L Z ). Same conventions as in 
Fig.H 



B. Energy J- (I) of the domain structure 

To simplify the expressions of the energy it is conve- 
nient to introduce the normalized lengths 
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where A e g = A 2 /L z is Pearl's penetration depth which, 
in the limit L z <C A, takes the place of A as the effective 
magnetic length scaled. Note that in these notations 
A c ff = A 2 . Furthermore we will consider the normalized 
energy JF = J-/(32ttMqL z ) and omit the bar henceforth. 
So the domain wall contribution is 7cw(i) = W/(32ttL) 
and the volume contribution, T vo i = Tu + -^sc is 
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where Q = (2k + l)ir/L and Q z = ^A- 2 + Q 2 . For 
a thick film, i.e. Q z = q z L z > 1 so that coth(Q^) = 
1, Faure and Buzdin's result 7 is recovered as expected. 
Since the limit I <C L z has been previously investigated 7 
we discuss below the limit I 3> L z , that is L ^$> 1, with 
more details. 



1. Regime A 3> / 

This parameter regime is realized in the normal phase 
where A = oo, as well as inside the superconducting 
phase, in the vicinity of the critical temperature T c . 
When A 3> I (i.e., when A 3> L), the sum in (ITU)) may be 
expanded in powers of L/A before summing. Hence for 
A> I > L z , 



while for I <C L z , 
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2. Regime I 3> L z and I ^> A c 



domain wall on a length scale A (compare for example 
the contour |B|/47rM = 0.4 in (a) and (b) of Fig. HJ. 
Note that the supercurrent is responsible for a kink that 
the lines of constant |B| show at the film surface (see 
Fig. H;b) and Fig. W^)). This is because V x B is dis- 
continuous at the interface when A~ 2 7^ 0. Flux lines also 
have a kink at the domain wall since the magnetization 
M is discontinuous there. 



This parameter regime may be realized in the super- 
conducting phase, for either thin films or large domains. 
When L>1 the sum in (I10[) can be approximated by an 
integral. Euler-MacLaurin's approximation yields 
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where the 1(A) is defined by Eq. (|A1|) in the Appendix. 
Then using asymptotic expression (IA2[) of 1(A) in the 
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FIG. 4. (color online). Normalized domain width ls/L z as a 
function of the normalized penetration depth X/L z for differ- 
ent values of the normalized effective domain wall thickness 
W = w/L z . 



limit A>1, one finds 
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and for A -C 1, asymptotic expression (|A3|) yields 
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III. DOMAIN WIDTH AT EQUILIBRIUM 
A. General results 

The equilibrium size Is is obtained by minimization of 
the total energy F(l) = J 7 dw(0 + -?Vol(0- Figures HHS] 
show Z5 as a function of the penetration depth A for dif- 
ferent values of the normalized wall thickness W . In the 
superconducting state A decreases with decreasing tem- 
perature from infinity at T = T c to a finite value at T = 0. 
As discussed below, just below T c the domain width al- 
ways decreases from the normal state value In realized 
for A = 00. One can distinguish two regimes. When 
the domains in the normal state are wide (In 3> L z ) the 
decrease is negligible, while when the domains are nar- 
row (In <C L z ) their width can drastically shrink down 
to Is ~ A. Then, in both regimes, Is reaches a min- 
imum value before diverging at a critical lower bound 
A c . This limit corresponds to the situation where the 
ferromagnetic induction is completely screened by the 
supercurrent^r— and there is no stray field that need to 
be accommodated by a DS. 

The A-dependence of the width is similar to the 
prediction for DS in superconductor-ferromagnet (S/F) 
bilayer a 12 i 13 except for two aspects. First the shrinkage 



FIG. 5. (color online). Normalized domain width ls/L z as a 
function of the normalized penetration depth X/L z for differ- 
ent values of the normalized effective domain wall thickness 
W. Is converges to the normal-state width In when A tends 
to infinity, and it diverges at the critical value A c . In and X c 
are increasing functions of W. 
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FIG. 6. (color online). Normalized domain width ls/L z as a 
function of the normalized penetration depth X/L z for differ- 
ent values of the normalized effective domain wall thickness 

w. 
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. Ex- 



factor in a S/F hybrid is limited and Is > y 
periments performed on heterostructures made of a ferro- 
magnetic garnet layer combined with a superconducting 
layer of Pb^ or Nbi£ have observed the shrinkage of the 
domain size due to superconductivity. The smallest ob- 
served shrinkage factor ~ 0.3 was substantially smaller 
than the theoretical prediction y2/3. This circumstance 
is a non-equilibrium effect and is related to the special 
dynamics of the vortex structure created in the supercon- 
ducting layer and coupled to ferromagnetic domains 15 . 
The second point is that although the DS can also be 
unstable in S/F bilayers when A decreases, it exists a 
range of parameters (when the thickness of the F-layer is 
approximately larger than half the domain width in the 
normal state) in which lg remains finite even in the limit 
A = 13 , in contrast with SFM where the DS is never 
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FIG. 7. (color online). Normalized domain width In/L z in 
the normal state as a function of the normalized effective do- 
main wall thickness w/L z . Dashed lines are the plots of the 
asymptotic formulas. 



stable for a vanishing A. 

Hereafter, complementary to the exact numerical solu- 
tion, we present the analytical results obtained in various 
asymptotic regimes. 



1. Normal state and vicinity of the transition 

In the limit L ^> 1, using expression (|lip valid for the 
vicinity of the superconducting transition where A 2 ^> L, 
one finds that the minimum of the energy F(L) is realized 
at 
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which yields the domain width 
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in the normal state (see Fig. [7]) . The domain width Ls 
just below the superconducting transition is then related 
to the normal-state width Ln by 
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In the other limit L <C 1 , the approximation (fT2|) yields 
the standard resul t 16 i 17 



LN ~\ 
(see Fig. |7|) and 
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in accordance with Ref. 7 . Hence the domain width al- 
ways decreases when the system enters the superconduct- 
ing phase. 



FIG. 8. (color online). Normalized critical penetration depth 
A c /L z as a function of the normalized effective domain wall 
thickness w/L z . Dashed lines are the plots of the asymptotic 
formulas. 
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FIG. 9. (color online). Normalized critical penetration depth 
A c /L z as a function of the normal domain width In/L z . 
Dashed lines are the plots of the asymptotic formulas. 



2. Stability of the periodic domain structure 

The periodic domain structure is always stable in the 
normal state where the energy minimum is obtained for 
a finite width L. This is not the case in the supercon- 
ducting state. In the limit L 3> A 2 and L> 1, the 
saddle-point equation for the energy is 



L 2 = 



2 A 4 



ttA - 27(A) - W/8 
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where the function 1(A) is defined by Eq. (|Aip in Ap- 
pendix. There is an energy minimum at a finite L only for 
ttA- 27(A) -W/8 > 0. Since (ttA- 27(A)) increases with 
A, this condition defines a lower bound A c below which 
the periodic structure is unstable. In the thought ex- 
periment where A decreases from infinity (normal state) 
to zero, the period Ls first decreases from Ln but then 
it increases back before diverging as (A — A c ) -1 / 2 (see 
Fig. 01). Note that the limit of vanishingly small A has 
been considered in Refi£ and the conclusion of the ab- 
sence of DS was made. When L^ ^> 1 in the normal 
state, the approximations of 7(A) yield 
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while when the normal-state width Ln ^ 1, the periodic 
structure exists down to 
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Figures [8] and [9] show that these analytical relations fit 
well the numerical results in the limit regimes. 

One can now understand the difference of shrinkage 
factor between a narrow domain and a wide domain. 
For a narrow domain with I and A < L z , the width 
in the superconducting state decreases with the penetra- 
tion depth as Is ~ A. The shrinkage goes on until the DS 
is destabilized by an excessive superconducting screen- 
ing. The minimum value that the width can have is then 
l™ lm r*j X c ~ l 2 N /L z , which is much smaller than Zjv- In 
the opposite limit when I and A> L 2 , the actual mag- 
netic length scale is A e ff, not A. And the domain width 
behaves as Is ~ A e ff when A e ff is smaller than Iff. Thus 
the minimum Z<? ln ~ A 2 /-L z ~ In, so that a wide domain 
does not significantly shrink while the penetration depth 
decreases to the instability value A c . 



B. The case of triplet superconductivity 

As we will see in a next section, the temperature de- 
pendence of the domain equilibrium size can provide a 
phenomenological way of distinguishing triplet from sin- 
glet pairing. Here, we describe the general formalism; nu- 
merical applications to some specific compounds will be 
shown later as examples, in the framework of Uranium- 
based superconductors. 



1. Modification of the wall energy due to the depletion of 
condensation energy 

In the case of a triplet superconductor with Cooper 
pairs fully spin polarized by the intrinsic magnetization 
of the ferromagnetic domains, a depletion of the order 
parameter occurs on a width ~ £ at both sides of the 
interface separating two domains of opposite magnetiza- 
tions. This inhomogeneity results in a loss of conden- 
sation energy ~ 2£ ccm d£,d y L z at a domain wall. At a 
first approximation the domain wall energy for T < T c is 
modified by replacing the effective domain wall width w 
by the temperature dependent w t defined by 



U0)/(L/4) 



w t 



(l+^(l -i) 3/2 ) 



where t = T/T c , and 
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is proportional to the ratio of the superconducting energy 
at T = to the magnetization energy. Here, we have 
assumed f (*) = f(0)/vT^1, A(i) = X(0)/y/l^t, and 
£cond(£) = £cond(0)(l — t) 2 . Using the asymptotic expres- 



sions relating w to the normal state width Iff, one can 
estimate u from experimental data in the limit Iff <C L z 
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FIG. 10. (color online). Schematic temperature dependence 
of the domain width Is for the singlet case (thin solid lines) 
and for the triplet case (thick solid lines) in (i) a film of thick- 
ness L z and (ii) in another one with a smaller thickness Lz/4. 
The dashed-line curves show the full A-dependence of Is in 
the case of a constant effective wall thickness (a) w = wt=i 
and (b) w — u>t=o- 
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and in the limit Iff 3> L z with 
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2. Behavior of the domain width below the transition 

A triplet order parameter yields a vanishing correction 
to the domain wall energy at the transition. The tem- 
perature dependence of the domain size is in the limit 
£jv< 1, 
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So, as in the singlet case, Is always decreases with de- 
creasing temperature in the vicinity of T c . However at 
a lower temperature and for a large parameter u), the 
behavior for the triplet case can be significantly differ- 
ent from the singlet case. Since for a given Lff there 



is an universal curve Ls(A) for the singlet pairing, it is 
possible to detect an experimental signature of triplet 
superconductivity as a deviation from this curve. This 
is schematically illustrated in Fig. QJJ] which shows the 
temperature dependence of L$ in the singlet and in the 
triplet case. When temperature is reduced from T c (on 
the right side where X(T) — > oo) to T = (on the left side 
where X(T) = A(0)), the Lg-curve in the triplet case in- 
terpolates between the curve (a) obtained with a constant 
w = w t= i and the curve (b) for w = w t=0 . The triplet- 
case curve follows the singlet-case curve with w = Wt=i at 
high temperature but at a lower temperature a deviation 
appears which is maximum at T = 0. 

In the expression (|2T)|) of L$ for the triplet pairing, the 
first-order term is proportional to L 2 while the second- 
order term is proportional to L~ 2 so one can expect that 
the singlet/triplet discrepancy can be amplified when re- 
ducing the thickness of the film. As shown in Fig. [TU] this 
discrepancy is enhanced when the ratios X(0)/L z and 
w t /L z become larger. For small enough thickness the 
difference is not only quantitative but is also qualitative 
as the triplet-pairing curve can increase with decreasing 
temperature while the singlet-pairing curve decreases. 



where the Lerch transcendent &l{z, s,a) 
Y,T=o zk /( a + k Y- K furthermore I > X 2 /L z 
pression (|31j) can be simplified as 



W, 



w„ 



1- 



2A 2 



ex- 



(32) 



Thus W sc increases with the width of the domain and, 
when the latter becomes large compared to the magnetic 
radius of a vortex, the work is simply equal to the cou- 
pling energy W m = 2L Z ® M . 

In the opposite regime A ^> I, the screening current is 
vanishingly small so its work is only equal to a fraction 
of Woo. For wide domains (I ^$> L z ), the dependence of 
W» r on the domain width is linear as 
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4GWoo I 

n 2 Anf 



(33) 



where Catalan's constant G ~ 0.916, while for narrow 
domains (I <C L z ), the dependence is quadratic as 
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IV. CONDITION FOR VORTEX APPEARANCE 



B. Interaction energy of a vortex-antivortex chain 



In this section we estimate the energy for the creation 
of a vortex-antivortex pair in the middle of adjacent do- 
mains. This additional energy is the sum of (i) the energy 
decrease when (anti)vortices are driven by screening cur- 
rents away from domain boundaries to domain centers 
and (ii) the concomitant increase of the interaction en- 
ergy between one vortex and its antivortex as they move 
apart from each other. For simplicity we neglect the 
vortex-core energy assuming the limit of large k = A/£. 



A. Energy decrease due to the screening current 



The force of the screening current acting on a vortex 9 
is f£ J-l [V x (b — 47rM)] x e z dz, where the quantum 
of flux $o = 2.07 x 10~ 15 T.m 2 . When the vortex moves 
from the domain wall to the domain center, the current 
then produces a work equal to 
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(30) 

where Woo = 2L z <f>oMo is the coupling energy of the 
flux carried by one vortex with the magnetization Mo of 
a ferromagnetic domain. 

The work has simple analytical expressions in limit 
cases. In the regime I ^> A, 
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The interaction energy of a chain composed by one vor- 
tex (antivortex) located at the center of every domain of 
positive (negative) magnetization is the magnetic energy 
required to separate by a distance I vortices and their 
paired antivortices created at domain walls. In the case 
A 3> I or A c ff 3> I, the energy per one vortesii has the 
asymptotic expression 
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We have assumed £ <§C I. In a large-K superconductor 
with £(0) ^C I, there is always a temperature regime close 
to T c in which A > I > f . 

In the opposite limit, for I 3> A e g S> L z the interaction 
energy per vortex is 
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while for L z ^> I ^> A, it is 
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C. Condition for the appearance 

From the above results we can derive the condition of 
the vortex stability in a superconducting domain struc- 
ture. The latter is assured when the total energy per vor- 
tex £ t t = £ vv — yV sc is negative. This condition defines a 
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critical value of magnetization M v above which vortices 
appear. One can distinguish four regimes in which M v 
has a simple analytical expression. In a film with narrow 
domains (i.e. I -C L z ), 



4ttM v 



$o , A 



4ttA 2 VC 



In - when I > A, (38) 



47rA4 = ^-|ln(- ) when / .V i :->,!>) 



and 



Note that Eq. (|38|) recovers the expected result that vor- 
tices appear when the magnetization AitMg exceeds the 
lower critical field of a non-magnetic superconductor de- 
fined as H* x = ^fx? hi/4). However, when the penetra- 
tion depth is larger than the domain width, the coupling 
of the vortex flux with the DS magnetization is not en- 
ergetically optimum so the critical magnetization Q39[) is 
higher than the mono-domain value H*± . 

I 



In the opposite limit of wide domains (i.e. I >• L z ), 
the critical magnetization is 



while 



47rM " = ^ ln (^f L ) i0r/ A "' ! - 



47r ^-T^ ln (f lh,,/ W 



(40) 



(41) 



In this regime of parameters, ln(A c ff/£) > In (A/£) so 
Eqs. (|40| and (|4T|) can be interpreted as an extension of 
results Eqs. ([38]) and ([39]) for the limit I > L z where A 
is replaced by A e ff. 



V. APPLICATION TO URANIUM-BASED 
SUPERCONDUCTING FERROMAGNETS 





2ttM 2 (J/m 3 ) 


m m 


A(0) H 


L z (mm) 


In W 


w (nm) 


£cond (J/m 3 ) 


U) 


UGe 2 


14200 


15 


1 


2 


4 


13.6 


40 


0.039 


UCoGe 


51 


15 


1.2 


0.15 


2 


45 


29 


2.36 


URhGe 


3070 


18 


0.9 


0.2 


20 


3450 


30 


6.5 xlO" 4 



TABLE I. Experimental values the magnetization energy 2-kMq, the coherence length £(0), the magnetic penetration depth A(0), 
the half-thickness L z and the normal-state domain width l N in the compounds UGez^— , TJCoGe 3 -^^^ and URhGe2»i£i22i24. 
The effective domain wall thickness w is obtained from the experimental value of Zjv and L z . For an estimate of the condensation 
energy, we used the BCS formula £ C ond ~ 0.166T C AC where AC is the specific heat jump at T c . The parameter ui is calculated 
with definition (1251). 



According to experimental data (see Table [IJ), the 
Uranium-based superconductors UGe 2 , UCoGe and 
URhGe are in the limit Ijf <C L z . The effective do- 
main wall thickness w can then be easily calculated in 
this limit from the experimental value of the domain 
width In in the normal state. We roughly estimate the 
zero-temperature condensation energy with the BCS for- 
mula £ C ond(0) ~ 0.166T C AC where AC is the volumic 
specific heat jump at T c . Among the three compounds 
UCoGe has the smallest magnetization resulting into the 
largest u ~ 1 (see Table QJ. It is then the most promising 
candidate where to look for possible observation of the 
non-monotonic temperature dependence of the domain 
width in the superconducting phase. Using its estimated 
parameters we plotted the temperature dependence of 
the domain width for the singlet and the triplet scenario 
in a film of half-thickness L z — 0.15 mm (Fig. ITT]) and 
L z = 3fim (Fig. [12). For L z = 0.15 mm (Fig. [IT) the 
domain width decreases below T c in both the singlet and 
the triplet case but there is a significant difference of 
size at T — between the two. By reducing the film 



I 

thickness the discrepancy is amplified. For instance with 
L z = 3fim (Fig. [T2) one obtains a qualitative difference: 
for the triplet case the width increases by 32% from its 
value in the normal state while for the singlet case it 
shrinks by 21%. 



VI. CONCLUSION 

We have done the complete analysis, within the Lon- 
don approximation, of the domain width Is at equilib- 
rium in a SFM film of arbitrary thickness 2L Z . We have 
shown that the ratio ls/L z follows a universal depen- 
dence on the normalized penetration depth X/L z and 
the normalized effective wall thickness w/L z (see Figs.0]- 
[5). In addition to the exact numerical dependence, we 
have derived analytical expressions of this relation in 
limit cases. In particular, we have recovered the previ- 
ously published results by Faure and Buzdin 7 which had 
been established for L 2 > A and Is, and we have com- 
plemented them with the analysis in the opposite limit 



1.5 



0.5 



w=45nm, X(0)= 1.2 |im 
L z =150 |im 



/ 



co=2.36 



*co=0 



0.2 



0.4 0.6 0.8 

T/T r 



FIG. 11. (color online). Temperature dependence of the do- 
main width Is (thick line) calculated with experimental pa- 
rameters for UCoGe in a film of half-thickness L z — 150 (jm. 
The thin line is ls{T) obtained with uj = 0, which corresponds 
to singlet superconductivity. 



0.38 



0.34 




0.26 



0.22 



FIG. 12. (color online). Temperature dependence of the do- 
main width Is (thick line) calculated with experimental pa- 
rameters for UCoGe in a film of half-thickness L z — 3 pm. 
The thin line is ls(T) obtained with uj = 0, which corresponds 
to singlet superconductivity. 



L z < A and Is- 

We have found that the domain width always decreases 
with temperature when the DS enters the superconduct- 
ing phase. The screening supercurrent induces a decrease 
of the domain width if the latter is larger than the pene- 
tration depth. However, the paramagnetic screening sup- 
presses the DS below a critical value A c of penetration 
depth, which means that the system may become mono- 
domain at a finite temperature if A(0) < A c . The domain 
shrinkage is relatively small when the domain width In in 
the normal state is much larger than the film thickness. 
Indeed, in this limit, the width decreases as A 2 /L z which 
is only a fraction smaller than In since A c < y ' L z In- 
With the other limit shape of domains, i.e. In ^C L z , 
the shrinkage can be important since then the domain 



width decreases as A while the critical penetration depth 
^c ~ ^Ar/^2 "C In- 

Furthermore, we have investigated the effect of the 
triplet pairing within an effective temperature-dependent 
renormalization of the domain wall energy. Actually, 
the DS magnetization alternatively suppresses the order- 
parameter components with opposite spin-projections. 
This inhomogeneity results in an additional positive con- 
tribution to the domain wall energy, that is absent for the 
singlet pairing. The supplementary term vanishes at the 
transition temperature but can be significant at a lower 
temperature. The variation of the domain width is then 
no more described by the universal dependence of ls/L z 
on X/L z obtained for the singlet case. This gives a phe- 
nomenological way of distinguishing triplet from singlet 
pairing in experiments. And, as we have shown, the dis- 
crepancy can be amplified by reducing the film thickness 
(see Figs. HU and H2). 

We have also established that vortices can appear when 
the magnetization exceeds a critical value 4nM v . In nar- 
row domains (i.e. I <C L z ), the latter is equal to the 
lower critical field H* ± of a non-magnetic superconduc- 
tor if A <C I, but in the opposite limit A 3> I, it is much 
larger than H^. For wide domains (i.e. I >■ L z ), we have 
found that the critical magnetization behaves similarly, 
with the difference that A is then replaced by Pearl's ef- 
fective penetration depth A e ff. 

The available experimental data (Table Q] and 
Refi.-3ii8-.24) fo r UGe 2 , URhGe and UCoGe show that 
in these compounds I <C L z and the conditions for the 
vortices appearance (section lTV]) are fulfilled. This means 
that the effects of domain shrinkage (considered in sec- 
tion lllll for the case where the vortices are absent) should 
be weakened. Another difficulty to experimentally ob- 
serve the evolution in the superconducting state may be 
related with the pinning of the domain wall and/or vortex 
pinning. At present there are no convincing experimental 
data on the change of the domain structure below T c . 

It has been recently observed in S/F bilayer o 14 ' 15 that 
the coupling between vortices and magnetic domains 
leads to a strong shrinkage of domains in the presence 
of the oscillating field used for equilibration of the do- 
main structure. It would be interesting to perform the 
similar experiments with superconducting ferromagnets 
to study these coupling effects. 

We thank K. Hasselbach, C. Paulsen, J. -P. Brison, and 
D. Aoki for helpful discussions. This work has been sup- 
ported by the French ANR project SINUS. 



Appendix A: The function 1(A) 

We define the function 1(A) as the integral 

u du 



1(A) 



o I 1 



< ; J-ucothu/^r 



(Al) 
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1. Limit A > 1 

The integral can be rewritten as 

T(A) - f~ Jtiu 

J V i4 J — J0(A- 2 +m 2 )(A- 2 +u 2 +m) 
„oo u 2 (l-VA~ 2 +u 2 coth \/A~' 2 +v?)du 

JO (A- 2 +« 2 )(A- 2 +« 2 +«)(A- 2 +« 2 +«VA- 2 +m 2 coth^/A^+u 1 ) ' 



SO 



by 



using 



VA~ 2 + u 2 coth VA~ 2 + u 2 



the approximations 

a 1 + (A" 2 + u 2 )/3 in 



the numerator and ss 1 in the denominator when u < 1, 



and VA -2 + u 2 



coth VA 2 + u 2 w u when w > 1, 



(iw. 



r/n ^ r» m <jm If 1 

I \ 1Y ) ~ J0(A- 2 +« 2 )(A- 2 +« 2 +«) 3 JO (A- 2 +« 2 +«) 2 



+ /" 



u (1 — u)du 



(A- 2 +u 2 )(A- 2 +u 2 +u)(A- 2 +2u 2 ) ■ 



Then for A > 1 



J(A)w — -21nA- — +ln2 
v y 2 12 



2. Limit A < 1 



(A2) 



In the limit A <C 1, coth VA 2 + u 2 w 1 then 

/(A) 



A z 



v / A ^ + « 2 ' 




(v^+« 2 +«) 


poo 


ii cfo 



'o Vl + v 2 (Vl + f 2 +v) 
(l-ln2)A 2 . 



(A3) 
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